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Chapter 1 



Introduction 



1.1 Why nonequilibrium Green functions ? 

In these lectures notes we discuss a method that is not very familiar to most quantum chemists, the 
nonequilibrim Green function method. However, this is likely to change mainly due to new devel- 
opments in nanoscience and the emerging field of 'molecular electronics' in which single molecules 
are used as conducting elements attached to electrodes. Description of such situations is a chal- 
lenge for the theorist as it requires a description fast time-dependent processes in strong external 
fields. The traditional quantum chemistry approaches that take into account electron correlations 
in stationary sytems can not deal with these situations. This has been our main reason for looking 
at the nonequilibrium Green function method, which has had important applications within solid 
state, nuclear and plasma physics. However, due to its general nature it can equally deal with 
molecular systems. Let us briefly describe its main features: 

• The method has as its main ingredient the Green function, which is a function of two space- 
time coordinates. From knowledge of this function one can calculate time-dependent expec- 
tation values such as currents and densities, electron addition and removal energies and the 
total energy of the system. 

• In the absence of external fields the nonequilibrium Green function method reduces to the 
equilibrium Green function method which has had imporatnt applications in quantum chem- 
istry. 

• Nonequilibrium Green functions can be applied to both extended and finite systems. 

• The nonequilibrium Green function can handle strong external fields nonperturbatively. The 
electron-electron interactions are taken into account by infinite summations. 

• The approximations within the nonequilibrium Green function method can be chosen such 
that macroscopic conservation laws as those of particle number, momentum and angular 
momentum are automatically satisfied 
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• Dissipative processes and memory effects in transport that occur due to electron-electron 
interactions and coupling of electronic to nuclear vibrations can be clearly diagrammatically 
analyzed 

1.2 References 

Many more things can be said about nonequilibrium Green functions. We therefore give a list of 
references for further reading that we found useful ourselves: 

Nonequilibrium Green functions 

• P.Danielewicz, Quantum Theory of Nonequilibrium Processes, Annals of Physics, 152, 239, 
(1984) 

• H.Haug and A.-P.Jauho, Quantum Kinetics in Transport and Optics of Semi-conductors, 
Springer- Verlag, Berlin (1998) 

• L.P.Kadanoff and G.Baym, Quantum Statistical Mechanics: Green's Function Methods in 
Equilibrium and Nonequilibrium Problems, Addison- Wesley, (1989) 

• M.Bonitz, Quantum Kinetic Theory, Teubner, Stuttgart-Leipzig, (1998) 

• J. Rammer and H.Smith, Quantum field-theoretical methods in transport theory of metals, 
Rev.Mod.Phys. 58, (1986) 

• G.Stefanucci and C.-O.Almbladh, Time- dependent partition-free approach in resonant tun- 
neling systems, Phys.Rev.B69, 195318 (2004) 

Equilibrium Green function theory 

• E.K.U. Gross, E.Runge and O.Heinonen, Many-Particle Theory, Adam-Hilger, Bristol (1991) 

• A. L. Fetter and J.D.Walecka, Quantum Theory of Many-Particle Systems , Dover, (2003) 

• L.Hedin and S.Lundqvist, Effects of electron- electron and electron-phonon interactions on the 
one-electron states of solids, Solid State Physics 23, 1969 

• J.Linderberg and Y.Ohrn, Propagators in Quantum Chemistry, Wiley-Interscience, (2004) 
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Chapter 2 



Second quantization 



2.1 The Schrodinger equation 

The goal of these lectures notes is to describe many-electron systems in general time-dependent 
external fields. We will restrict ourselves to a nonrelativistic description and therefore we restict 
ourselves to a discussion of the time-dependent Schrodinger equation (TDSE). The TDSE of n- 
particles has the form 

(id t -H(t))*(x 1 ,...,x n ,t) = (2.1) 

where Xj = riat denotes a space-spin variable of particle i. Now the Hamiltonian is invariant under 
interchange of two particles. This implies that the eigenstates of the Hamiltonian can be chosen to 
transform according to an irreducible representation of the permutation group. Experience teaches 
us that only the one-dimensional representations, i.e. the completely symmetric or completely anti- 
symmetric representations seem to occur in nature. If this assumption is made, a famous theorem 
of Pauli states that the symmetric representation necessarily describes particles with integer spin 
whereas the anti-symmetric representation necessarily describes particle with half-integer spin. The 
proof is based on the properties of the representations of the Lorentz group and is therefore an 
essentially relativistic result. It would go to far to describe this in detail and we simply use the 
result. Since electrons are spin-half particles they are described by anti- symmetric wave functions: 

. . Xi . . . Xj . . .) = -*(. . . Xj . . . Xi . . .) (2.2) 

In the following section we will introduce a formalism that will automatically take this symmetry 
into account within the operators. This formalism is commonly knwon as second quantization and 
is used in most advanced many-body approaches. Knowledge of this formalism is therefore essential 
to understand a great number of research papers. 
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2.2 Fock space and field operators 
2.2.1 Definitions 

In this section we will define a space, the so-called Fock space, of quantum states that consist of 
linear combinations of states with different number of particles. There are several reasons for doing 
this. 

• This will provide us with a natural framework to discuss processes that change the number 
of particles, such as in ionization. 

• It will allow us to consider temperature dependent systems and calculate expectation values 
within the grand canonical ensemble in which the particle number is not fixed. 

• By dealing with creation and annihilation operators in Fock space the anti-symmetry prop- 
erties of the wavefunctions are automatically built into the anti-commutation relations of the 
operators and we never have to deal with the wavefunctions themselves. 

Let {\^j,N}} be a complete set of states in a iV-particle Hilbert space. Then Fock space is defined 
as the set of linear combinations of the form 

oo oo 

\F) = a \0} + J2J2^\^j,N) (2.3) 

N=l j=\ 

The zero-particle Hilbert space is one-dimensional and has only one basis function |0) which can 
be indentified with the number 1. Let 

oo oo 
N=l j=l 

be another element of Fock space with respect to some other basis {|$j,jv)} then the inner product 
is defined as 

oo oo 

(F\G) ee al(i, + £ E a f* (*jM*i,N) (2.5) 

JV=1 i,j=l 

where (^jvI^jv) is the inner product in iV-particle Hilbert space. One can check that definition 
Eq.(2.5) satisfies all the usual requirements of an inner product. We also see that according to our 
definition the Fock overlap between states containing different number of particles ia zero, i.e. if 
\F) = |*i,jv) and \G) = \^j,m) then (F\G) = when N / M. Since we will only consider fermions 
we further assume that all the iV-particle Hilbert spaces consist of anti-symmetric functions. 
We are now ready to define the annihilation and creation operators that map wavefunctions from 
an n-particle Hilbert space to, respectively a n — 1 and a n + 1-particle Hilbert space. We define 
the action of the annihilation operator as 

(V)(x)*„)(xi,...,X n _i) EE v / " , I'n(xi,...,X n _i,x) (2.6) 



6 



This operator maps an n-particle state to a n — 1 particle state. If ^ n is anti-symmetric it is not 
difficult to see that the resulting n — I-particle state is again anti-symmetric. Corresponding to 
this annihilation operator V>(x) we define the adjoint operator ^'(x) which maps a n-partcle state 
to a n + I-particle state. This operator is called the creation operator. The explicit action of this 
operator on a n-particle state \l/ n is 

(^t( x )* n )( Xl , • • • , x n+1 ) = ^=J= ^(-iy +1 <5(x - Xj )* n ( Xl , . . . ±j . . . x n+1 ) (2.7) 

where Xj denotes that this argument is lacking. For example (xi, &2, X3, X4) = (xi,X3,X4). Let us 
check that ffl is indeed the adjoint of V>. We have 

<*„_i|^(x)|* n ) = v 7 " J dxi...x n _i**_ 1 (xi...x n _i)* n (xi...x n _i,x) (2.8) 

On the other hand we also have 

n 

(v^Wl^n-i)* = (-lrTt-l)^ 1 /dx 1 ...x n * n (x 1 ...x n ) 

x <5(x - x J -)$*_ 1 (xi . . . Xj . . . x n ) 

1 n f 

= -^(-l)"- 1 V(-1) J+1 / dxi...x n * n (xi...x J -_i,x,x J -+i...x n ) 

Vn 3=1 J 

X**_!(xi.. 

• x i • • • x nJ 

1 n r 

= (-!)"-! ^(-l)^(-l)^- / dx 1 ...x n * n (x 1 ...x j _ 1 ,x j+1 ...x n ,x) 

j=l 

X $n-l( X l • • -Xj_i,Xj + i . . .X n ) 
= 7^5Z / dx l--- dx n-l*n-l( x l--" X n-l)*n(xi...X n _i,x) 

= v^y'dxi...dx n _i$*_ 1 (xi...x n _i)* n (xi...x ri _i,x) (2.9) 

From Eqns.(2.8) and (2.9) we see that 

<* n _i|^(x)|* n ) = (* n |^(x)|$ n _i)* (2.10) 
From this result it is simple to prove that for any Fock space states F and G we have 

W(x)|G) = <G|^t( x )|F>* (2.11) 
This proves that $'(x) is the adjoint of $(x). 



Exercise 

Prove relation (2.11) from Eq.(2.10) 
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We did not check sofar that ■0^(x)|* T j) is indeed an anti-symmetric n + 1-particle state if * n is 
antisymmetric. However, this is readily checked as well. Let us interchange elements p and q with 
p > q. Then we first write 

^(x)|tf„> = 

+ 
+ 

The first term in this equation is certainly anti-symmetric if we interchange p and q. Let us therefore 
consider the last two terms. Interchanging x p and x 9 for these terms gives 

(-l) p ^ n (xi . . .x ? _i,Xp,x g+ i . ..Xp. . .x n+ i) + (-l) p * n (xi ...-k q ,. . .Xp_i,x 9 ,Xp+i . . .x n+ i) 
= (-1)? (-l)^ +1 *„( Xl . . . X, . . . x n+1 ) + (-l)«(-l)^ +1 * n ( Xl . . . Xp . . . x n+1 ) 

= -(-l)9* n (xi...x,...x n+ i) - (-l)f* n (xi...x p ...x n+ i) (2.13) 

which indeed yields an overall minus sign. We thus see that the creation operator maps an anti- 
symmetric n-particle state to an anti-symmetric n + 1-particle state. 

2.2.2 Anti-commutation relations 

We subsequently prove a basic anti-commutation relation satisfied by the field operators. We 
calculate 

^ t (x)^(x')^„ = V> t (x)v / ^1'„(xi,...,x n _i,x / ) 

n 

= (-1)"- 1 ^(-ly+^x - x,-)*„(xi ... x, ... x n , x') 

n 

= E(- 1 )' +1<5 ( X " ^)(-l) n - j *n(x!, . . • , X^x', Xj . . . X n ) 

3=1 

n 

= ^«5(x-Xj)*n(xi,...,x J -_i,x',Xj...x„) (2.14) 

J'=l 



(-1)' 

(-1)" 

(-1)" 
y/n + 1 



n+l 



V; £ (-l^x-x^xi,...^- 
v n + 1 ■ n 7i ^ 

(-l) p ^( Xl ...x p ...x n+1 ) 
(-l) 9 *(xi...x 9 ...x n+ i) 



X n+ lj 



(2.12) 
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We further have 

71+1 



j>(x!)ft(x)* n = ^)_^(-l)"^(-l)^ 1 ,J(x-x J 0*„(xi,...x J -...,x n+ i) 

Vn + l ^ 

n 

= ("I)" £(-l) J ' +1 <S(x " ^)*n(xi X„,x') 
J'=l 

+ (-^"(-lJ^X-xO^X!,...,^) 

n 

= (-l)"X;(-l) J ' +1 (-l) n -^(x - X,)* n ( Xl . . .X^X'^ . . .X n ) 

J'=l 

+ 5(X - X / )*n(xi, . . . ,X n ) 

n 

= - ^ <5(x - Xj)*„(xi . . . Xj_i, x', Xj+i . . . x n ) 

+ 5(x-x')*n(xi,...,x n ) (2.15) 
From Eqns.(2.14) and (2.15) we find 

^(x')^(x) + Vit( x )^( x ')j ^ n = § (y L - x ')*„ (2.16) 
for any wave function We thus obtain the basic anti-commutation relation for the field operators 

{^t( x ),^(x')} = *(x-x') (2.17) 
where anti-commutation brackets are defined as 

= AB + BA (2.18) 
Further anti-commutation between the field operators are readily derived. We have 

■0(x)'0(x / )*n = \/Wn - l* n (xi . . .X n _ 2 ,X,x') 



= -yfn\Jn- I* n (xi . . .x n _ 2 ,x',x) 

= -V>(x')^(x)* n (2.19) 

for any ^> n and therefore we have 

{^(x),^(x')}=0 (2.20) 
By taking the adjoint of this relation we then immediately also have 

{^ ( x)^t (x ')} =0 (2 . 21) 

2.3 The Hamiltonian in second quantization 



With the relation derived in the previous section it is not difficult to express the Hamiltonian in 
terms of field operators. Any one-body multiplicative or differential operator O can be written as 

6 = J^o(xi) = J dx$t(x)o(x)^(x) (2.22) 
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Proof. Using Eq.(2.14) we have 

j dx^( x ) ( x )^( x )|* n ) = j dxo(x / )V' t (x)^(x / )| x;=x |^) 

= ^2 J (ixo(x / )(5(x-x j )^ n (xi...x i _i,x / ,x j+ i...x n ) 

3 

= J>(x,)|*„> (2.23) 

j 

For example for the second quantized expressions for the kinetic energy and a time-dependent 
external scalar potential are given by 

f = Y.-\^ = -\j rfx^(x)V 2 ^(x) 

i 

V(t) = Y^v(xit) = J dx^ t (x)^(x)t;(x,t) (2.24) 

i 

Another example, is, for instance, the density operator 

Hy) = Y, S & ~ x *) = J dx<5 ^ " x)^(x)^(x) = ^(y)^(y) (2.25) 

i 

We can further consider two-body operators. Let us consider the expression 

n(, x i, ■ • ■ > x n— 1) y> X J 

rt— 1 

= y^(-l)^ 2 V;t( x ) J2(-l) j+1 S(y - x,)* n ( Xl , x n _ 1; y , x) 

3=1 
n—1 

= v^(-l) n ^(x) X;(-l) J ' +1 (-l) n_1 "^(y - X,)*n(x!, . . . ,X n _!, X) 
3 

= v^W^W-Xn-i) (2.26) 
where we defined the wavefunction 

4?(xi...x„-i) = [<J(y-xi) + ... + < y(y-x n )]* n (xi...x n _i,x) (2.27) 

Now since 

n 

$^i ) (xi...x jfc ...x„) = ^ <5(y - Xj)$ n (xi . . .x fe . . . ,x n ,x) 

= (-l)"- fe *(y-x J -)*„(xi...x fc _i,x,x fc+ i...,x n ) (2.28) 

it follows from Eq.(2.26) and Eq.(2.28) that 

n 

= (-I)"" 1 ^(-l) fe+1 5(x - Xfc)*^(xi ...x; ...x„! 

fc 

n 

= Yl ^(x-x fe )5(y-x j )* n (x 1 ,...,x n ) (2.29) 
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So we derived the expression 

n 

^ t (x)^t(y)^(y)^( x ) = ^ <5(x - Xjfc )<5(y - Xj ) (2.30) 

We see that for the two-particle interaction we can write 

Y n If n 

^ = -J2 w (^^k) = ^dxdyw(x,y)^25(x-x k )5(y-x j ) 

= l -J d«fy«;(x,y)^t( x )^t( y )^( y )^( x ) (2.31) 

Therefore a general Hamiltonian of the form 

n \ n 

H(t) = J2H*j,t) + 2 I>( x i> x *) ( 2 - 32 ) 

can in terms of the field operators be written as 

H(t) = J dxft(x)h(x,t)fa) + iy dxdy W (x,y)^t( x )^t( y )^( y )^( x ) (2 .33) 

This is the main result of this section. Let further derive some other often used representation 
of the Hamiltonian. Let </>i( x ) be a complete set of orthonormal orbitals in a one-particle Hilbert 
space, i.e. we have 



Sij = J dx^*(x)<^(x) (2.34) 

oo 

<5(x-y) = X>*(x)^(y) (2.35) 



i=l 

Then we can define the annihilation and creation operators with respect to this basis as 



Oi = J dx99*(x)^(x) (2.36) 

4 = j rfx^(x)Vit(x) (2.37) 

From the properties of the complete set of orbitals and the anti-commutation relations of the field 
operators we then immediately find the expressions 

oo 

V>(x) = ^2<pi(y)ai (2.38) 

i=l 

oo 

Vit(x) = j>?(x)aj (2.39) 
i=l 

Sij = {a\,a 3 } (2.40) 

= {a i ,a j } = {a\,a]} (2.41) 
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In terms of the operators hi and a] the Hamiltonian the attains the following form 



H(t) = M*) a\aj + V vki a\a)a k ai (2.42) 



i,j i,j,k,l 



where we defined the coefficients 



hij(t) = 




/ d*dy < ^*(x)^(y) u; ( x 'y)^fe(y)^( x ) 



(2.43) 



(2.44) 



Exercise 

Derive the form of the Hamiltonian in Eq.(2.42). 



2.3.1 Some useful relations 

Finally we will derive some useful relations that play a role when we will discuss the equations of 
motion for the Green function later on. These relations are derived with straightforward use of the 
anti-commutation relations of the field operators. Let us evaluate the commutator 



$(x),^(y)^( Z )] 



V'(x)^ t (y)V'(z) - $t(y)$(z)$(x) 
"<5(x - y) - ^(y)^(x)] ^(z) - ^(y)^(z)^) 




= (S(x-y)^(z) 



(2.45) 



Let O now be a one-body operator. Then we immediately see from Eq.(2.45) that 




z=y 



(2.46) 
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This was our first useful relation. Let us continue with the commutator of a two-body operator: 

$(x),^(y)^t( Z )^( Z )^(y)] = $(x)^(y)^t( Z )^( Z )^( y ) _ ^t(y)^t( Z )^( Z )^(y)^( X ) 



<J(x - y) - ^(y)V(x)j^ T (z)V'(z)^(y) - ^ ' {yW r (z)V(z)^(y)V'(x) 
= «5(x-y)Vit( z )^( z )^(y) -^t (y) ^ (x) ^t (z) ^ (z)] +^t (z) ^ (z) ^ (x) ^ (y) 

-^(y)^(z)^(z)^(y)^(x) 
= <5(x - y)^(z)^(z)^(y) - 5(x - z)^( y )^( z )^( y ) _ ^t (y) ^t (z) ^ (z) |^ (x) ^ (y) j 
= 5(x-y)^( z )^( z )^( y ) + ( j( x _ z )^t( y )^( y )^( z ) (2.47) 

With this relation we see immediately that the commutator of the field operator ^(x) with the 
two-particle interaction W is given by 

$(x),W] = I J dydzw(y,z)$(x)JHy)$Hz)j>(z)$(y)] 

= \j ^y dzw (y' z )(^( x -y)^ t ( z )^( z )^(y) + ^( x - z )V ,t (y)V'(y)V'( z )) 

= J dzu-(x,z)^ t (z)^(z)^(y) (2.48) 
This is our second useful relation. 



Exercise 

Derive analogously the following relations 

[V>t( x ),V>t (y) ^ (z)] = -<5(x-z)Vit( y ) (2.49) 
[^(x), ^t( y )^t( z )^( z )^( y )] = _j( x _ y)^t( y )^t( z )^( z ) _ «5( X _ z )^( z )^t( y) ^( yI 2.50) 

and consequently the following expressions for the commutator with a one-body operator O and 
two-particle interaction W: 

[^(x),6] = -o(x)^t( x ) (2.51) 
[^(x),^] = -J dzw(x,z)^(x)V> t (z)^(z) (2.52) 



We now collect our results and calculate the commutator of the field operators with the Hamiltonian. 
We the relations derived above we obtain: 

[^(x).fl-(t)] = h(x.,t)$(x) + Jdyw(x.,y)$\y)4(y)fa) (2.53) 

[^(x),ff(t)] = -h(x,t)4,\yi)- J dyw(x,y)$\x)ft(y)4,(y) (2.54) 



13 



Chapter 3 

The time contour 



3.1 The evolution operator 
3.1.1 Definition 

Let |*(t)) be a solution to the TDSE then the time evolution operator U(t,t') is defined by the 
relation 

m)) = uM\W)) (3-i) 

i.e. it maps a wave function at time t' to a wave function at time t. It obviously satisfies the 
relation U(t,t) = 1. If we differentiate Eq.(3.1) with respect to t and use the TDSE we obtain 

id t U(t,t')\V(t')) = H(t)U(t,t')\V(t')) (3.2) 

Since | is arbitrary we find that 

id t U(t,t') = H(t)U(t,t') (3.3) 

On the other hand if we differentiate Eq.(3.1) with respect to t' we obtain 

= (id t ,U(t,t'))Mt')) + U(t,t')H(t')Mt')) (3.4) 

This must be again true for any \^f(t')) and therefore we have 

idfU(t, t') = -U(t, t')H(t') (3.5) 

We thus see that the evolution operator satisfies the relations 

id t U{t,t') = H(t)U(t,t') (3.6) 
id t >U{t,t') = -U(t,t')H(t') (3.7) 
U(t,t) = 1 (3.8) 
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These relations completely define the properties of the evolution operator. If we integrate Eq.(3.6) 
from if to t where t > tf we have 



U(t,t')-U(t',t') = [ dt 1 d h U(t 1 ,t') = -i [ dhH(h)U(h,t') (3.9) 
Jt' Jt' 

and we obtain the expression 

U(t,t') = l-i f dhH(h)U(h,t') (3.10) 
Jt' 



if 

This expression can be iterated to give 

ft ftl rtn-i 



^ ft ftl ft n -l 

U(t, t') = 1 + J^(-i) n / dh / di 2 ... \ di n H(h) ■ ■ ■ H{i n ) 
n= 1 

oo / '\n ft ft ft 

= l + Y.^T- / d *i / dt 2 ... di n T[H(h)...H(i n )\ (3.11) 
^ n! Jf Jt' Jt' 



where we in the last step could make all integrals run from t' to t by introducing the time-ordered 
product 

T[H(h) . . . H(i n )\ = ^2 9(i P{1 ) - t P(2 )) . . . 0(ip( n -i) - ip(n))H(tp(i)) ■ ■ ■ H{ip(n)) (3-12) 
p 

where P runs over all permutations of the numbers 1 . . . n. The latter expression is known as the 
time-ordered product which orders the operators at the latest times to the left. Similarly we can 
integrate Eq.(3.7) from t to tf when tf > t. This gives 

U(t,t')-U(t,t)= [ di 1 d h U(t,i 1 ) = i [ dhU{t,h)H{h) (3.13) 
Jt Jt 

and therefore we have 

U(t,t') = l + i J dhU(t,h)H(h) (3.14) 
Iteration of this equation now gives 

f*t (*t~l f t n \ 

U(t, t') = l + ^i n dh / dt 2 .-. di n H{i n ) . . . H(h) 



n=l 

00 jn rt ft ft 

= 1 + V -. / dh di 2 ... di n f[H(h) . . . H(i n )} (3.15) 
n ~r 1 n - Jt' Jt' Jt' 

where now we defined anti-chronological time-ordering by 

f[H(h) . . . H{i n )\ = 0(i P{ i) ~ i P (2)) ■ ■ ■ 0(i P (n-i) ~ i P( n))H(ip {n) ) . . . H(i P{1) ) (3.16) 
p 

i.e the latest times are ordered to the right. We can this formally write 

U(t,t') = Texp(-i / dTH{r)) when tf < t (3.17) 
Jt' 

U(t,t') = fexp(ij^ dTH{r)) when t 1 > t (3.18) 
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3.1.2 Evolution of ensembles 



In the remainder of these lectures we will discuss the evolution of systems that are originally 
in thermodynamic equilibrium. In that case statistical mechanics tells us that the expectation 
values of operators are given as traces of appropriate ensembles. Let us therefore recall some basic 
definitions. The trace of an operator A is defined as 

Tr {i} = (3-19) 

i 

where is a complete set in Fock space. The trace is independent of chosen set. This is easily 

seen by insertion of a complete set: 

= ^(^l^^lil^^^^-lil^) (3.20) 

i,j 3 

If we now consider a system at temperature T then its equilibrium expectation value of operator 
O in the grand canonical ensemble is given by 

where (3 = 1/ksT and fee is the Boltzmann constant and where present energy eigenstates of 
the system. This expression can equivalently be written as 

(6) = Tr jp6} (3.22) 

where we defined 

(3.23) 



When we apply an external field to a system initially in thermodynamic equilibrium, it will subse- 
quently evolve in time. It is this time-dependennce that we want to study. Our Hamiltonian will 
have the general form. 

H(t) = j dx${*)h{r,t)i>{yL)+ 1 -J J d Xl dx 2 ^t( Xl )^t( X2 ) _l_ ^( X2 )^( Xl ) ; (3.24) 
where 

h(r, t)= l - HV + A(r, t)] 2 + v(r, t) - M . (3.25) 

In the latter equation we introduced external potential v(r, t) and vector potential A(r, t) which 
are switched on at t = to, i.e for t < to our Hamiltonian will be time-independent. Note that we 
also included the chemical potential in the definition of h. This means that we can simplify the 
expression of the statistical operator p to 

P=- W (3-26) 

Tr e-P H " 
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Figure 3.1: Keldysh contour. The forward and backward parts of the contour are on the real axis 
but are plotted slightly off this axis to display the two branches more clearly. 

where we denote by Hq the Hamiltonian for times t < to- After we switch on an external field the 
initial equilibrium ensemble will evolve in time and the new expectation value becomes 

^I ^My ^lAW} (3.27) 

where 

d H (t) = U(t ,t)d(t)U(t,t ) (3.28) 

is the operator in the Heisenberg picture. Our task is now to evaluate such expectation values. 
First we will analyze the expression Eq.(3.27) a bit further. 



3.2 Contour ordering 
3.2.1 Definition of the contour 

The operator e~* H ° can now be regarded as an evolution operator in imaginary time, i.e. 

U(to - i0,to) = e~ pA ° (3.29) 

if we define H(t) to be equal to Ho on the contour running straight from to to to — i@ m the complex 
time plane. We can therefore rewrite our expression for the expectation value as 

Tr \u(t -iP,t )U(t ,t)dU(t,t )\ 

W)) = 1 F 1 1 (3-30) 

Tr {U(t -i(3,t )\ 

If we read the time arguments of the evolution operators in the numerator of this expression from 
left to right we may say that the system evolves from to along the real time axis to t after which 
the operator O acts. Then the system evolves back along the real axis from time t to to and finally 
parallel to the imaginary axis from to to to — i/3. A corresponding contour, origianlly introduced 
by Keldysh [1], is displayed in Fig. 3.1. With this observation we rewrite Eq.(3.30) as 
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Tr { T c [exp ( -i L dtH{t) ) 6 (<)] } 

W)) = 1 — ( : y 1 (3.31) 

Tr [T c [exp(-if c dtH(t))}} 



where we define 



T c [exp (-i [ dtH(t))6(t)} = V / dh . . . dinTcM^Hih) . . ■ H(i n )} (3.32) 

Jc ±^ nl J c 



n=0 

and we further defined the contour-ordered product 



Tc[M{t\) . . . A n (t n )\ — 0(*P(1),*P(2)) • • • ^(*P(ti-1))*P(ti))^P(1)(*P(1)) ■ ■ ■ ^P{n){tp(n)) (3.33) 



p 



where 6{t\,t2) are contour step functions 9 generalized to arguments on the contour [2], 

f 1 if t\ is later than ty on the contour 
0(tl ' tl,) = \o otherwise (3 ' 34) 

For instance, time t\ in Fig. 3.1 is later than time t2 on the contour. All the time-integrals in 
Eq.(3.32) are taken along the contour. 



3.2.2 Functional derivatives and the time-dependence of expectation values 



This observation motivates us to define the following action functional 

S = ilnTr {t/(t - */Mo)} > ( 3 - 35 ) 
where we define the evolution operator on the contour as 

U(t,t') = T c exp(-i / dtH{t)). (3.36) 
Jt> 

Let us now see how this functional can be used as a generating functional by making variations 
with respect to parameters in the Hamiltonian. To do this one needs to consider changes in U. 
When we make a perturbation 5V(t) in the Hamiltonian we have using Eqs.(3.6), (3.7) and (3.8) : 

id t SU(t,t') = SV(t)U(t,t') + H(t)SU(t,t') (3.37) 
id t >SU(t,t') = -U(t,t')SV(t') - 5U{t,t')H(t') (3.38) 
SU(t,t) = (3.39) 

The solution to this equation is given by 

SU(t,t') = -i [ dTU(t,T)SV(T)U(T,t') (3.40) 
from which variations in the action can be calculated. 
Exercise 

Check that Eq.(3.40) is a solution to Eqs.(3.37) and Eq.(3.38) satisfying boundary condition 
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Eq.(3.39). 

For instance, if we choose 



5V(t) = J dx5v(xt)h(x) (3.41) 

where n(x) = i/>t(x)i/i(x) is the density operator we have 
8U(to-ip,t ) 



6v(xt) 

From this equation we have 

SS 



= -iU{t - ip, t)h(x)U(t, t ) = -iU(to - i/3, t )h H (xt) (3.42) 



Tr {u(t -i(3,t)h(yL)U(t,t )} 

Sv W Tr {U(to-ip,to)} 

Similarly we can also calculate the change in expectation values. 



<n(x,t)> (3.43) 



-(0(h))- ' J 



Sv(xL 2 t 2 ) 6v(xL 2 t 2 ) Tr {u(t -ip,to)} 

Tr {[/(to - %P, ti)6(ti)L>(ti, t 2 )h(x 2 )U(t 2 , t )} 

Tr {#(t -i/Mo)} 
Tr {[/(to - t 2 )h(yi 2 )U(t 2 , ii)d(ti)#(*i, t )} 
Tr {[/(t to)} 



= -iO(h,t 2 ) 

-iO(t 2 ,h) 



Tr {[/(to-i/3,ti)0(ti)[/(ti,to))Tr {[/(to - »/3, t 2 )n(x 2 )[/(t 2 , t )j 
+j L J L J. (3.44) 

Tr {[/(t to)} Tr [u{t - ip,t )} 

where the last term follows from differentiation of the denominator. The equation (3.44) can be 
rewritten as 

= -^(<i ) *2)(Oi/^i)nff(x 2 t 2 ))-^(t2,ti)(nff(x 2 t 2 )dH(ti))+i(0(ti))(n if (x 2 t 2 )) (3.45) 

di>(x 2 t 2 ) 

If we define the fluctuation operator 

AO^(xt) = 6 H (xt) - <On(xt)) (3.46) 

this can be rewritten as 

= -^(ti,t 2 )(A6 ff (t 1 )An ff (x 2 t 2 ))-^(t 2 ,t 1 )(An J/ (x 2 t 2 )Ad if (t 1 )) 

du(x 2 t 2 J 

= -i(T c [Ad H (h)Ah H (x 2 t 2 )}) (3.47) 
Similarly we can calculate derivatives of time-ordered products. We have 

5{TC foteh) t2)]) = - i ^[ i ^)^^)^^3)]) + ^^(^3))(Tc[(i(tl)B(t 2 )]) 

= -i(T c [i(ti)5(t 2 )An // (x 3 t 3 )]) (3.48) 
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Exercise 

Prove relation (3.48) 

We finally consider the time-dependence of expectation values and consider the time-derivative of 
an expectation value (0(t)). We then can calculate using Eqs.(3.6) and Eq.(3.7) : 

idt6 H (t) = idt U(to,t)6(t)U(t, t ) = -U(to,t)H(t)6(t)U(t, t ) 
+U(t , t)6(t)H(t)U(t, t ) + iU(to, t)dtO(t)U(t, to) 

= u(t ,t){[d(t),H(t)] + id t d(t)}u(t,t ) 

= [d H (t),H H (t)]+i(d t d(t)) H (3.49) 

where in the last step we used that for two operators A and B we have 

U(t ,t)A(t)B(t)U(t,t Q ) = U(t ,t)A(t)U(t,t )U(t ,t)B(t)U(t,to) 

= A H (t)B H (t) (3.50) 

When we use expression (3.49) in the definition of the expectation value (3.30) we obtain 

id t {6{t)) = ([6 H (t),H H (t)}) + i((d t 6(t)) H ) (3.51) 

3.2.3 Calculations with contour-ordered functions 



In this section we will derive some relations that we will use later on. For simplicity we will 
introduce the notation i = x±ti for the space-time variables. In the previous section we saw that 
differentiation of expectation values naturally led to the consideration of contour-ordered operators 
of the form: 

<r c [A(l)B(2)]) = e(ti,* 2 )<A(l)5(2)) + (9(* 2 ,t 1 )<B(2)A(l)) (3.52) 
We therefore consider general functions of the form 

a(t, t') = a s {t) S(t, t') + 0(t, t') a>(t, t') + 6(t', t) a<(t, t') (3.53) 

where the delta function on the contour is defined as dt9(t,t'). A function of this form is said to 
belong to Keldysh space. When we consider the equation of motion of the Green function later we 
have to calculate contour integrals of products of such functions, i.e. function of the form 

c(M') = / dia(t,t)b(i,t') (3.54) 
Jc 

where both a and b belong to Keldysh space. We will see that such a function also belongs to 
Keldysh space. Let us work out this integral. We have : 

c(M') = a{t,t')b 5 {t') + a 5 {t)b{t,t') 

+ / di \e(t, i)6(t, i')a> b> + 9(t, t)6(t', t)a> b< + 0(t, t)6(t, t')a < b> + 6(t, t)6(t', t)a> 6 > 1(3.55) 
Jc 1 J 
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We can now consider several cases let t > t' on the contour and be real. Then we have 

c>(M0 = a > (t,t')b s (t') + a 5 (t)b > (t,t') 

+ f dia > {t,t)b > (t,t')+ [ dta > (t,t)b < (t,t')+ [° " dta < (t,i)b > (i,t') (3.56) 
Jt' Jto Jt 

We denote 

a < (t,t -ir) = a\t,r) (3.57) 
a > (t -iT,t) = a r (r,i) (3.58) 

This notation is quite suggestive when one reads the symbols ] and [ from left to right. For instance 
] has a horizontal segment followed by a vertical one; correspondingly a) has a first argument that 
is real (and thus lies on the horizontal axis) and a second argument that is imaginary (and lies on 
the vertical axis). With the notations of Eq.(3.57) and Eq.(3.58) we can rewrite expression (3.56) 
as 



c>(M') = a > (t,t')b 5 (t') + a 5 (t)b > (t,t')+ [ dt^^t) -a < (t,i)]b > (t,t') 

Jt 

rt' l-P 

- / dta > (t,t)[6 > (t,t / ) -b < (i,t')]-i dTa\t,T)b^{T,t') (3.59) 
Jt Jo 

When we define the retarded and advanced functions in physical time as 

a R (t, t') = a 5 (t)5(t - t') + 9{t - t')[a>(t, t') - a < (t, t')} (3.60) 
a A (t,t') = a 5 (t)5(t-t') -6(t' - t)[a>(t, t') -a < (t,t')} (3.61) 

then Eq.(3.59) can be written shortly as 

■(t,t')= dia R (t 1 t)b > (t,t')+ dia>(t,t)b A (i,t) -i (ira 1 (^,T)fe ^ (r,^ , ) (3.62) 
Jt Jt Jo 



c>< 



This expression can be even more simplified when we introduce the notation 

poo 

a-b = / dta(t)b(t) (3.63) 
Jt 

a*b = -i dra(T)b(T) (3.64) 
Jo 

Then we simply find 

c > = a R ■ b> + a> ■ b A + aU (3.65) 

From Eq.(3.56) we can also derive an expression when t > t' on the contour and when t is on the 
imaginary part of the contour t = to — it. Equation (3.55) then becomes 

c r( r , t') = c>(to - ir, t') = a > (t - ir, t')b 5 (t') + a 5 (t - ir)5 > (t - *r, t') 

rto-ir rt' 

+ / dta > (t -iT 1 t)b > (t,t')+ dta > (t - ir,i)b < {t,t') 

Jt 1 Jto 

+ / dta < (t,t)b > (t,t / ) (3.66) 

J to — IT 



21 



This expression can be rewritten as 

rt' 

c^(t,0 = a^T,t')b\t') + a\t -iT)b^T,t')- dta^r,t)[b > (i,l/)-b < (i,t')] 

Jto 

/•to— ir rto—if3 

+ / dia > {to-iT,t)tf{i,t') + dia < (t - iT,t)tf(f,t') 

Jto JtQ — ir 

poo p(3 

= / dia^T,i)b A (i,t')-i dTa M (T,f)tf(T,t') (3.67) 

Jtn JO 



where, for convenience, we defined the Matsubara function a M as the function a restricted to the 
imaginary part of the contour 

« M (n, r 2 ) = i a s (h)S(n - r 2 ) + 0(n - r 2 )a > (i 1 ,i 2 ) + 0(r 2 - r 1 )a < (t 2 , h) (3.68) 

where in this expression t\ = to — iT\ and i 2 = to — ir 2 . Equation (3.67) can in condensed notation 
be rewritten as 

c r = J-6 A + a M *6r (3.69) 

Finally, if we let t and t' be both on the imaginary part of the contour then one readily sees that 
the real part of the contour does not contribute in Eq.(3.54) and we immediately find that 

M M , .M 



c M = a M *b M (3.70) 



Exercise 

Derive the relations 



c< = a R -b< +a< •6 A + a 1 *6 r (3.71) 
c l = a R -b l +aUb M (3.72) 
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Chapter 4 

The Green function 



4.1 Definition 



Before we give the definition of the Green function we first extend our definition of time-ordering 
slightly. Sofar we only considered operators at a particular time that consist of products of an 
even number of field operators, such as the Hamiltonian and operators that represent observables. 
Such operators we call non-fermionic (they are, however, not bosonic in general since they need 
not satisfy bosonic commutation relations). However, in the following we will also consider time- 
ordering of fermionic field operators themselves. To treat this case the definition of time-ordering 
is extended to 

Tc[M{h) ■ ■■A n (t n )\ = ^ (-l) Fp 6'(tp( 1 ),tp(2)) . ■ ■0(?P(n-l),*P(ri))4p(l)(*P(l)) • • • ^P{n){tp(n)) 

(4-1) 

where Fp is the number of times that a fermionic operator gets interchanged with another fermionic 
operator in the permutation P. Let us give an example. Let O be a non-fermionic operator that 
consist of products of an even number of field operators. Then according to our definition we have 

T C [^(1)^ (2 )0(3)] = 

9(h,t 2 ) 6(t 2 , t 3 ) ^(1)^(2)6(3) + 6(h,t 3 ) 6(t 3 , t 2 ) ^(1)6(3)^(2) 
+ 6(t 3 , h) e(h,t 2 ) 6(3)^(1)^(2) - 9(t 2 , h) 9(h,t 3 ) $,(2)^(1)6(3) 
- 0(t 2 , t 3 ) 6(t 3 , h) $,(2)6(3)^(1) - 9{h, t 2 ) 9(t 2 , h) 6(3)^(2)^(1) (4.2) 

We see that the last three terms gained a minus sign as ^(2) and ^p-(l) were interchanged in 
these terms. With this definition of time-ordering we now define the one-particle Green's function 
G as 

1 Tr {U(t - i0,t o )Tc fe/(l)V^(2)l ) 
1 Tr [U(t -i(3,t )\ 

= -i(T c [^(l)^(2)]), (4.3) 
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This can be written in the form 



G(xit - i0, 2) 



G(l,2) = e(t 1 ,t 2 )G > (l,2) + 0(t 2 ,t 1 )G < (l,2). (4.4) 

where we introduced the greater and lesser Green functions G > and G < according to 

G>(1,2) = -i(fa(l)ft H (2)) (4.5) 
G<(1,2) = i(ft(2)^(l)) (4.6) 

Let us now derive the boundary conditions that G satisfies. If we consider the Green function at 
t\ = to — i(3 and use the cyclic property of the trace we find 

1 Tr {^( Xl )^ - i/3, t 2 )^ t (x 2 )f7(t 2 ,t )} 
* Tr [u(t -ifft )} 

x Tr (?7(to-i/3,t2)^ t (x 2 )C>(t 2 ,to)V'(xi)) 

= t L fz s i - = -G(x 1 t ,2). (4.7) 

1 Tr [U(t -ip,t )\ 

The Green function defined in Eq. (4.3) therefore obeys the boundary condition G(xiio,2) = 
— G(xiio — £2)- The property G(l,x 2 io) = — G(l,x 2 io — iff) for the other argument is likewise 
easily verified. These boundary conditions are sometimes referred to as the Kubo-Martin-Schwinger 
conditions [3, 4, 5]. Similar boundary conditions are satisfied by the usual equilibrium temperature 
Green function which, in fact, is obtained for the special case where the time arguments are located 
on the contour along the imaginary axis to to to — where the Hamiltonian is time-independent. 
Analogously to the one-particle Green function one can further define the n-particle Green function 

G„(l, . . . , n, 1', . . . , ri) = (-i) n {T c $ H (l) . . . ^(n)^(l') • • • (4.8) 

Also these Green function satisfy Kubo-Martin-Schwinger boundary conditions. One can derive a 
set of so-called hierarchy equations of motion that relate the n-particle Green function to the nil 
-particle Green functions. 



Exercise 

Derive the condition 

G(l,x 2 to) = -G(l,x 2 to-i/?) (4.9) 



4.2 Physical content 
4.2.1 Expectation values 



The Green function G is a basic ingredient in a diagrammatic perturbation expansion and it directly 
gives us the expectation values of one-particle operators. A one-body operator O can be expressed 
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as 

(O(t)) = /dxo(x't)(V>t(x,t)^(x't)) = -i [ dxo(x't)G < (xt,x't) (4.10) 
7 x '= x J x '= x 

For instance, if we denote by 1 + the limit to t± from above on the contour, the density is quite 
simply 

<n(l)) = -zG(l,l+) (4.11) 



and the current density is 



<j(l)> = ~i 



Vi 

2i 



2i 



+ A(1) 



g(i,i'; 



(4.12) 



i'=i+ 



The calculated observables will obviously depend on what approximation scheme we use to obtain 
G. It is therefore important that these approximations are such that the calculated observables 
satisfy the macroscopic conservation laws, like e.g. the continuity equation, dt{h) = — V • (j). Such 
approximations are called conserving approximations and will be discussed later on. 



4.2.2 Removal and addition energies 

We will now derive a useful relation for the Green function that clearly displays its physical content. 
This relation is known as the Lehmann representation. We first consider the Green function on the 
imaginary part of the contour. In that case we have t = —it where r runs from to (3 (we put the 
time to to zero without loss of generality) . In that case we have 

G(xi, -rri,x 2 , -in) = G M (xm, x 2 r 2 ) 
= 6>(ti - r 2 )G > (xi, -in,x 2 , -in) + 6>(r 2 - n)G < (xi, -in,x 2 , -in) (4.13) 

Let us first analyze G K . We have 

G < (xi, -in,x 2 , -in) = i(V' t (x 2 , -ir 2 )^(xi, -in)} 

= i^ J](^|e-^ o e^ 0T2 ^ t (x 2 )e-^ ,T2 e^ 0T1 ^(xi)e-^ 0T1 |^) 

i 

= ^e-^e^^ 2 -^)<^|^ t (x 2 )e-^ ,T " 2 |^ i )<^ i |e^»^^(xi)|^ i ) 

Z — ^ 

h3 

= 4E e_/3£<e(f;i " E ' )(T "" Tl) (*^ t ( X2 )i^)^i^( xi )i^) ( 4 - 14 ) 

where we defined Z = Tr {U(to — i/3,to)}. We can do a similar calculation for G > 

G > (xi,-in,x 2 ,-ir 2 ) = -z(?/>(xi, -iri)^ t (x 2 , -in)) 
= -i\ V(^|e- /3i?0 e^ ,Tl V;(xi)e-^ oTl e^ 0T2 ^ t (x 2 )e-^ ,T2 |^) 

Zj — ' 

% 

= -i^^e-^e^( ri - T2 )(^|^(x 1 )e- i?0Tl |^)(^| e ^ r2 ^ t (x 2 )|^) 

- -4E e "^ le(£l " Bj)(ri " T2) ^^( x i)i^)^i^( x 2)i^) ( 4 - 15 ) 
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Let us check the anti-periodicity conditions. We have for < r 2 < (3 : 

G( Xl ,0,x 2 ,-zr 2 ) = G < (x 1 ,0,x 2 ,-ir 2 ) (4.16) 
G(xi,-z/?,x 2 ,-ir 2 ) = G > (xi,-i/3,x 2 ,-ir 2 ) (4.17) 

and thereforewe must have 

G < (x 1 ,0,x 2 ,-ir 2 ) = -G > (x 1 ,-i/3,x 2 ,-ir 2 ) (4.18) 

From the explicit expressions Eq.(4.15) and Eq.(4.14) we see that this relation is indeed satisfied. 
When we consider finite electronic systems we can often take the zero temperature limit {(3 — > 00). 
In that case we can choose the chemical potential such that Eq < and Ei > 0. If the ground state 
then has N particles we find for G> the following expressions 

G<(xi, -iri,x 2 , -ir 2 ) = i £ e^-^-^^ (^ \^ (y: 2 )\N - l,j)(N - 1, j|^(xi)|* > 

3 

G > (x 1 ,-ir 1 ,x 2 ,-ir 2 ) = -i^e^-^+ 1 ^( ri - r2 )(^ |^(xi)|iV + l,i)(A^ + 1, jI^^x,)!^) 

3 

where \N ± 1, j) denote iV ± 1 -particle eigenstates of the system. The calculation above could, of 
course, also easily have been carried out in real time on the real axis provided we do not switch-on 
any time-dependent external fields. In that case we have (for simplicity in the zero-temperature 
limit) 

j 

G > (xi,ti,X2,t 2 ) = -i^e^»-^+ 1 ^(* 1 -' 2 )(^o|^(xi)|^ + l,i)(iV + l,i|^ t (x 2 )|^o) 

3 

With the definitions of the so-called Feynman-Dyson aplitudes 

#(xi) = <iV- l,i|^( Xl )|* > (4.19) 
/,( Xl ) = (* |^(xi)|iV+l,j) (4.20) 

These expressions simplify to 

G < (x 1 ,x 2 ;t 1 - t 2 ) = i^eWw-^-^to-VgjMgfai) (4.21) 

3 

G > (x 1 ,x 2 ,t 1 -t 2 ) = -iY,e i{ENfi - EN+1 - j)(tl - t2) fjixi)f;(x2) (4.22) 

j 

We can now for the equilibrium situation define the spectral functions 

x 2 ; u>) = i f drG^(xi, x 2 ; r)e^ T (4.23) 

We find the explicit expressions 

A<(xi,x 2 ;w) = -^^(xi^CxaJ^w-^o-^-ij)) (4.24) 
3 

A>(xi,x 2 ;w) = ^/ J (x 1 )/;(x 2 )^ + (^o-i?iV+i J )) (4-25) 

i 

We thus see that Fourier transforming G< gives a function that is peaked at the ionization energies 
of the system, whereas Fourier transforming G > gives a function that is peaked at the addition 
energies or affinities. For instance, photo-electron spectra can directly calculated from knowledge 
of A<. 
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4.2.3 Excitation energies 



The time-propagation of the Green function gives also the expectation values of the time-dependent 
one-body observables. A subsequent Fourier transformation then gives us also directly the excita- 
tion energies of the system. Let us illustrate this with the density operator. The time-dependent 
density is directly obtained from the Green function from Eq.(4.11). If we expand the density in 
powers of the external field we obtain 

(Mi,) = M.) + / * + \ j |^«(^(») + ■ • ■ 

= no(xi) + (-i) J d2 {T c [Ah H (l)Ah H (2)})v(2) 

+ { ~2~ j ( T cl A nH(l)Ah H (2)Ah H (3)})v(2)v(3) + ... 

(4.26) 

where we used Eqs.(3.47) and Eq.(3.48). If we undo the contoour integrations we obtain the linear 
and higher order density response functions. Let us for instance consider the linear term. We write 

x (l, 2) = -i(T c [An H (l)An H (2)]} = 9(t 1 ,t 2 ) X > (l,2) + 0(t 2 , ti)x < (l, 2) (4.27) 

where 

X > (1,2) = -i(Ah H (l)Ah H (2)} (4.28) 

X > (1,2) = -i(An H (2)Ah H (l)} (4.29) 



Then we have 



<M1)> " no(xi) = f d2 X (l, 2)v(2) = T d2 X > (1, 2)«(2) + / ° X <(1, 2)«(2) 

= / d2(x > (l,2)- X < (l,2)M2)= / XR {l,2)v{2) (4.30) 
J to J to 

where we used that external perturbing field is zero on the imaginary part of the contour and we 
define the retarded response function as 

X r(1, 2) = 9(h - i 2 )(x > (l, 2) - x < (l, 2)) = -i9(h - t 2 )([h H (l),h H (2)}) (4.31) 

and where we used that the commutators of two operators is the same as the commutators of the 
corresponding fluctuation operators. Now the operators x > have a similar Lehmann representation 
as the Green function. For instance in the zero temperature limit we have 

X < (xi,ti,x 2 ,t 2 ) = -z^e^°-^)( t2 - tl )(*o|n(x 2 )|^,i)(iV,j|n(x 1 )|*o) (4.32) 

j 

This expression is very similar to that of G < . The main difference is that instead of removal 
energies this expression contains particle- number conserving excitations. Fourier transformation 
of this quantity gives us the excitation spectrum of the system. Similar Lehmann representations 
exist for the higher order response functions that will play a role when we propagate the Green 
function in stronger extrenal fields. 
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Chapter 5 



The equations of motion 



5.1 The self-energy 

We now study the equation of motion for the Green function. Using the definition of operators in 
the Heisenberg picture, and the Hamiltonian as given in Eq. (3.24), the equation of motion for the 
field operators is 

= h(l)j> H (l) + Jd2w(l, 2)^(2)^(2)^(1) (5.1) 

id tl ^ H (i) = [^ H (i),H H (h)] 

= -Ml)^(l)- ^2^(1,2)^(1)^(2)^(2) (5.2) 

where i«(l,2) = <5(<i, t2)/|i"i — r 2 \ is the Coulomb interaction. The notation S(t\,t 2 ) again indicates 
that the time-arguments are on the contour. In deriving these expression we used Eqs.(2.53) and 
(2.54) together with Eq.(3.49). Let us now take the time derivative of the Green function: 

id tl G{\,2) = id h e(h,t2)G > (l,2)+id ti e(t 2 ,t 1 )G < (l,2) 
+9(t 1 ,t 2 )id tl G > (l, 2) + 0(t 2 , ti)i^ tl G < (l, 2) 
= i5(t l ,t 2 )[G > (l,2)-G < (l,2)} 

-ie(ti,t 2 )<i^ 1 ^H(l)ft(2)) + ie(t 2 ,t 1 ){4>\ I (2)id t J H (l)) (5.3) 

For the first term in this expression we can use the anti-commutation relations at equal times 

iSituh)^ (1,2) - G<(1,2)} = afa.taX^xitiJ^xa,*!)}) = 5(1,2) (5.4) 

where 5(1,2) = <5(xi —x. 2 )5(ti,t 2 ). Using Eqs.(5.1) in the last two terms of Eq.(5.3) we then obtain 

iQ tl G(l,2) = (y(l,2) + /i(l)G(l,2) (5.5) 

- / / d3«;(l,3)[e(ti,t 2 )<ft(3)^(3)^H(l)ft(2))-«(t2,ti)<ft( 2 )ft( 3 )^( 3 )^H(l)) 
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Now in the last line of this equation the pair of operators $ H (3)^(3) = n#(3) is actually evaluated 
at at £3 = ti due to the presence of an equal-time delta function in w(l, 3). However, this pair of 
operators is always situated to the left of V>h(1) and this ordering can thus also be obtained from 



lim(r c [^(l)^(2)n H (3)]) 
0— >u 



t 3 =tl+<5 

t aw au. +.\/„/.t 



= 0(ti,t2)(n H (3)V'H(l)^(2))-e(t 2 ,*i)<^T r (2)riH(3)^H(l)) (5.6) 

*3=*1 

where t\ + 8 means that the limit is taken from above on the contour. This equation also follows 
directly from Eq.(4.2) by taking 0(3) = n#(3). If we then replace io(l,3) by w(l + ,3) where 
1+ = xi,£i + 5 we can write Eq.(5.5) as 

(i9t 1 -/i(l))G(l,2) = 5(1,2) -t y d2,w{l + ^){T c [ipH{l)^ H {2)h H {2,)}) (5.7) 

Using the definition of the two-particle Green function (4.8) we can also write 

(T c [4> H (l)ft H (2)n H (3)]) = G 2 (l,3,3+,2) (5.8) 
We can therefore write the equation of motion as 

[id tl - h(l)]G(l,2) = 5(1,2) -i J d3u>(l+,3)G 2 (l,3,3+,2). (5.9) 

Similarly one has the adjoint equation 

[-id t2 -h(2)]G(l,2) = 5(1,2) -i J d3^(2+,3)G 2 (l,3,3 + ,2). (5.10) 



Exercise 

Derive Eq.(5.10) 



The problem is that the equation of motion for G depends on the two-particle Green function G 2 . 
Instead of propagating the equation with some approximate form of G 2 , we introduce the electronic 
self energy E, such that the term —1G2W is replaced with EG. We define the self-energy E and its 
adjoint X by the equations 

y"d2E(l,2)G(2,l') = -i J d2w(l,2)G 2 (l, 2, 2+,l') (5.11) 

J d2G(l,2)t(2,l') = -i J d2w(l',2)G 2 (l,2,2 + ,l') (5.12) 

It is not difficult to derive that E = E for initial equilibrium conditions. To do this we first write 
Eqs.(5.1) and (5.2) in the following form: 

(id tl - h(l))4> H (l) = i H (l) (5.13) 

H3 tl -Ml))4(l) = *Jr(l) (5-14) 

where we defined the operators 

i H (l) = J d2w(l,2)4,\ I (2)$ H (2)$ H (l) (5.15) 

V H {1) = J d2w(l, 2)^(1)^(2)^(2) (5.16) 
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These operators are readily seen to be each others adjoints. In terms of these operators the equations 
of motion of the Green function become 

(id tl -h(l))G(l,2) = 5(1,2) -i{T c [i H (l)ft H (2)}) (5.17) 
(-id t2 -h(2))G(l,2) = 5{l,2)-i(T c $ H {l)V H {2)\) (5.18) 

where from the derivation it is clear that the operators in and must be regarded as fermionic 
in the time-ordered product. From these equations we see immediately that 

j d3E(l,3)G(3,2) = -i{T c [i H {l)^ H {2)}) (5.19) 

J d3G(l,3)E(3,2) = -i{T c $ H {l)%{2)]) (5.20) 

From the equation of motion of the Green function we see that for the left hand sides of Eqs.(5.19) 
and (5.20) we can write 

{-id t2 - h(2)) J d3X(l,3)G(3,2) = E(l, 2) + J d3 dA E(l, 3)G(3, 4)E(4, 2) (5.21) 

(id tl - h(l)) J d3G(l,3)t(3,2) = E(l,2) + J d3d4E(l,3)G(3,4)S(4,2) (5.22) 

whereas for the right hand sides we obtain 

-i(-id t2 -h(2))(Tc[i H (l)^ H m) 
= -d t2 9(h, t 2 )(i H {l)$ H {2)) + d t2 e(t 2 , ^(^(2)^(1)) - i(T c [^(l)^(2)]) 

= <5(t 1 ,t 2 )({^(l),^(2)})-z(T c [^(l)^(2)]) (5.23) 

and 

-i{id tl -h{l))(T c $ H {l)V H {2)}) 
= d tl 9{t u t 2 )$ H {l)V H {2)) -^^2,^)^(2)^(1)) -i(T c [i H (l)iU 2 )}) 
= <5(ti,t 2 )({^(l),i{ r (2)}>-i<Tc[iH(l)it(2)]) (5-24) 
We have therefore obtained the relations 

E(l,2) + y d3d4E(l,3)G(3,4)E(4,2) 

= <5(t 1 ,t 2 )({^(l),^(2)})-i(T c [^(l)4(2)]) (5.25) 
E(l,2) + y d3cZ4E(l,3)G(3,4)E(4,2) 

= 5(*i,t 2 )<{^H(l),i{ r (2)})-i(ro[?H(l)it(2)]> (5-26) 

Now it is readily seen by computation that for the equal time parts on the left hand sides of 
Eqs.(5.25) and (5.26) we have 

<{i ff (xiti), ft( x 2*i)}> = <{^ H (xiti),ik(x 2 ti)}> (5.27) 
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and therefore Eqs.(5.25) and (5.26) immediately imply that 

£(1,2) = £(1,2) (5.28) 
We further see from Eq.(5.25) that £ has the following structure 

£(1,2) = £ 5 (1,2) +0(t 1 ,t 2 )£>(l,2) + 0(t 2 ,ti)E<(l,2) (5.29) 
where the equal-time singular part of £ is given by 

£ 5 (1, 2) = 6(h, t 2 )({iH(xiti), ft(x 2 ti)}) (5.30) 

This term is reaily worked out in more explicit form using the equal time anti-commutation relations 
of the field operators 

|iij(xiti),^(x 2 ti)| 

= J dx 3 w(xi, x 3 ) ^^(x 3 ti)^ // (x 3 ti)^H(xiti)^(x 2 ti) + ^(x 2 ti)^(x 3 ti)Vi i j(x 3 ti)'0 jH -(xiti)^ 
= <5(xi-x 2 ) y"dx 3 ^(xi,x 3 )^(x 3 ti)^(x 3 ti) - ^(xi,x 3 )^(x 2 ii)^(xiti) (5.31) 
We therefore obtain 

£ 5 (1,2) = 5(1,2) y'dx^X^XgKf^Xgt!)^^ 

= -z<5(l,2) J dx 3 w(x 1 ,x 3 )G < (x 3 t 1 ,x 3 t 1 ) + i5(t 1 ,t 2 )w(x 1 ,x 2 )G < (x 1 t 1 ,x 2 t 1 ) (5.32) 

This expression will later be denoted as the Hartree-Fock form of the self-energy. 

In this section we showed that £ = £ and therefore the equations of motion can be written as 

[id tl -h(l)]G(l,2) = 5(1,2) + y*d3£(l,3)G(3,2). (5.33) 

[-id t2 -h(2)]G(l,2) = 5(1,2) + y*(23G(l,3)£(3,2). (5.34) 

We see later that if we consider more general initial conditions that £ is no longer equal to its adjoint 
£. The next task is find more explicit forms of the self-energy. The self-energy is a functional of 
the one-particle Green function, and as a consequence, Eqs. (5.33) and (5.34) constitute a set of 
equations that should be solved to self-consistency once the functional dependence of £ on G is 
known. This is the topic of the next section. 

5.2 Derivation of self-consistent equations 

To generate a set of self-consistent equations we first note that from Eq.(3.48) it follows that 

= "'^3) < Tc ^ h(1) ^ (2)1) = -(Tc[^H(l)^ H (2)h H (m + <^(3)><r c [^(l)^(2)]> 

(5.35) 
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and we can therefore write 

-t<r c [^(l)ft(2)n H (3)) = + (n H (S)) G(l, 2) (5.36) 

With this expression we can write the equations of motion Eqs.(5.9) and (5.10) as [6] 
[id tl -h(l)]G(l,l') = 8(1,1') 

+ i J d2w(l + ,2) 5G ^ ) +G(l,l') J d2w(l,2)(h H (2)) (5.37) 
[- i d t ,-h(l')]G(l,l') = 5(1,1') 

+ i J d2w(l' + ,2) 5G ^ ] +0(1,10 J d2w(l',2)(h H (2)) (5.38) 

As a next step we will derive an expression for 6G(l,l')/5v(2). To do this we differentiate the 
equations of motion Eq.(5.33) and Eq.(5.34) with respect to v and obtain the equations: 

- ft < 1 »»%f-/ rf3E < 1 ' 3 » W =i(1 ' 2)0(u ' ) + / d3 W GM (5 ' 39) 

H^-W)^|p-/<i3^ifE(3,l') = *(l',2)G(l,l') + / dSCd.Sl^Hp (5.40) 

From the equation of motion of the Green function we can see that a general solution to these 
equations is 

^ikp = G(l,2)G(2,lO + J d3d4G(l,3)^|^G(4,lO + C(l,l / ,2) (5.41) 
where C is a solution to the homogeneous equations 

(id tl - h(l))C(l, 1', 2) - J d3 E(l, 3)C(3, 1', 2) = (5.42) 

(-z5 t / - /i(lO)C(l, 1', 2)- J d3 C(l, 3, 2)E(3, lO = (5.43) 



Exercise 

Check that the expression in Eq.(5.41) is a solution to Eqs.(5.39) and (5.40). 



To determine the solution uniquely we must use the boundary conditions. Since both terms on the 
right hand side of Eq.(5.35) satisfy the Kubo-Martin-Schwinger boundary conditions see that 

SG^to-ip,!') <SG(xi,t ,lO 



5v(2) 5v(2) 
5G(l,id 1 ,t -if3) _ 5G{l,^ v t ) 
5v(2) Sv(2) 



(5.44) 
(5.45) 



It is readily seen from the boundary conditions on the Green function that the first two terms in 
Eq.(5.41) satisfy the boundary conditions Eq.(5.44) and (5.45). Therefore also the function C must 
satisfy these boundary conditions. This uniquely fixes C = as the solution to Eqns.(5.42) and 
(5.43). We thus have 



5G(1,1') 
5v(2) 



G(1,2)G(2,1') + J d 3d4G(l,3)^|^G(4,lO (5.46) 
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From the structure of this equation is then natural to define the following so-called vertex function 
T as 

r(12;3) = 5(l,2)5(2,3) + ^p (5.47) 

such that we have 

^j^p = j d3 d4 G(l, 3)G(4, l')r(34; 2) (5.48) 
If we insert this equation into Eq.(5.37) we obtain 
[id tl -h(l)]G(l,l')= 5(1,1') 
+ i J d2d3d4G(l,3)u;(l + ,2)r(34;2)G(4,l / ) + G(l,l') J d2w(l,2){h H (2)) 

= <5(l,l') + yd4E(l,4)G(4,l') (5.49) 
We therefore see that we can write £ as 

£(1,2) = z J d3d4G(l,3)w(l+,4)r(32;4) -i<5(l,2) J d3 w(l, 3)G(3, 3+) (5.50) 

where we used that (n#(3)) = — iG(3,3 + ). This equation can used iteratively to generate expres- 
sions for £ in terms of the Green function. This becomes more clear when we insert into Eq.(5.50) 
the explicit form of the vertex (5.47) : 

£(1,2) = iG(l,2)u>(l+,2) -iS(l,2) J d3w(l,3)G(3,3 + ) 

+ i j d3d4G(l,3)«;(l + ,4)^lp (5.51) 

The first line in this expression is defined as the Hartree-Fock self-energy 

£ HF [G,w](l,2) = iG(l,2)w(l + ,2)-iS(l,2) J d3w(l,3)G(3,3 + ) (5.52) 

The first iteration of Eq.(5.51) gives 

£(1,2) = £ HF (l,2)+i J d3d4G(l,3)^(l+,4) ^ ( ( 4 3 ) ' 2) 

+ z 2 | d3d4G(l,3)«;(l + ,4)^ j d5d6 G(3, 5)u;(3+ 6)^|^ (5.53) 

We can now evaluate (5£ HF /5v using Eq.(5.46): 

<5£ HF (3,2) . SG(3,2) f , ^G(5,5+) 

^(4) = ^^ir^ 3 ' 2) - i<5(3 ' 2) y d5 " (3 ' 5) ^r 

= zG(3,4)G(4,2)w(3+,2) — i<5(3, 2) ^ d5 u>(3, 5)G(5, 4)G(4, 5+) 

- i<5(3,2) / d5 rf6 rf7 w(3, 5)G(5, 6) G(7, 5+) 

J ov[A) 

+ iw(3 + ,2) j dh d6 d7 G(3, 6) G(7, 2) (5.54) 
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s = 




Figure 5.1: Diagrammatic expression of E^ 2 ) 

When we insert this expression back into Eq.(5.53) we obtain 

£(1,2) = E^(l, 2) -i 2 j d4d6d7 G(l, 2)w(l + , 4)w(2, 5)G(5, 6) G(7, 5+) 

+ i 2 / d3d4d6d7 G(l, 3)w(l, 4 + )w(3, 2+)G(3, 6) G(7, 2) 

J ov{4) 

+ i 2 [ d3d4G(l,3)w(l + ,4)t^ I d5d6 G(3, 5)W3+ 6)^%^ (5-55) 
J ov[4) J dv(Q) 

where we defined E< 2 ) as 

S( 2 )[G,^](1,2) = S HF (l,2)+i 2 y d3d4G(l,3)w(l + ,4)G(3,4)G(4,2)w(3 + ,2) 

- i 2 J d4d5G(l,2)u>(l + ,4)u;(2,5)G(5,4)G(4,5 + ) (5.56) 

We see that £( 2 ) for a given G and as functional of w is of second order in w. This term can be 
represented diagrammatically. To do this we associate with every Green function G(l, 2) a directed 
line running from 2 to 1 and with every interaction line io(l,2) a wiggly line connecting point 1 
and 2. This is illustrated in Fig. 5.1 for the expression given in Eq.(5.56). The remaining terms in 
Eq.(5.55) that contain derivatives of E are of higher order in w. For instance a further iteration of 
Eq.(5.55) produces new terms that are of third order in w. The main message is that we have found 
a systematic way of expressing the self-energy E as a functional of G and w. The whole iterative 
set of equations be written compactly by writing a coupled equation for the self-energy E and the 
vertex V. If we regard E as a functional of the Green function Eq.(5.47) implies immediately that 

= <S(1,2)<S(1,3) + J d4dM6d7 | G(4, 6)G(7, 5)r(67; 3) (5.57) 

We have therefore obtained the following set of self-consistent equations 

= [id tl -h(l)]G(l, 10-5(1,10 - J d2T,{l,2)G{2,\') (5.58) 

= [-ify -h(l') G(l,l')-*(1, l')-Jd2 G(1,2)E(2, l') (5.59) 

E(l,2) = i J d3d4G(l,3)w(l + ,4)T(32;4) -i5(l,2) J d3w(l,3)G(3,3 + ) (5.60) 

r(12;3) = <5(1,2)(5(1,3)+ /d4d5d6d7^4vG(4,6)G(7,5)r(67;3) (5.61) 

7 dG(4,5) 
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These equations can now be iterated to obtain self-consistent equations for the Green function. 
For instance, if we take the simplest approximation for the vertex, namely T(12; 3) = 6(1, 2)<5(1, 3), 
and insert it into Eq.(5.60) we obtain the Hartree-Fock approximation to the self-energy E^^. By 
inserting this expression into Eq.(5.61) we obtain a new approximation for the vertex from which 
one can obtain a new self-energy. As one can readily convince oneself the time-local part (i.e. 
proportional to , £2) ) of E is only given by T, HF , i.e. E 5 = E HF . It is therefore convenient to 
single out the Hartree-Fock part of the self-energy, such that [2] 

E(l,2) = E HF (l,2) + (9(*i,t 2 )E > (l,2) + e(t 2 ,*i)E < (l,2) (5.62) 

From expression (5.52) we see that we can write 

E HF (1, 2) = E HF (x 1 t 1 , x 2 ti)<$(*i, t 2 ) (5.63) 

where 

E HF (xi,ti,x 2 ii) = •iG < (xiii,x 2 ii)u;(xi,x 2 ) - i5(xi - x 2 ) J dx 3 w(xi, x 3 )G < (x 3 t 1 , x 3 ii) (5.64) 

Some specific forms for the self-energy will be discussed later. Let us first discuss the form of the 
equations of motion that one has to solve in any practical application. 



5.3 The Kadanoff-Baym equations 



Using Eq.(5.62) we transform the contour integration to obtain: 



(id tl -h(l))G^(l,2) = 1^(1,2) (5.65) 

(-id t2 -h(2))G$(l,2) = if (1,2) (5.66) 

(id tl -h(l))G\l,2) = /I (1,2) (5.67) 

(-id t2 -h(2))tf(l,2) = 1^(1,2) (5.68) 

(-d T1 -h(l))G M (l,2) = i5(n-r 2 )+I M (l,2) (5.69) 

(d T2 -h(2))G M (l,2) = ^(n-T 2 ) + / M (l,2) (5.70) 

where we defined the collision terms 

if = (EG)^ = E R -G^ + E^-G A + EUG r (5.71) 

if = (GE)^ = G R ■ E^ + G^ • E A + ★ E I" (5.72) 

/I = (EG) 1 = S R • G 1 + £l * G M (5.73) 

l\ = (GE)T = G r -E A + G M *E^ (5.74) 



The Eqs. (5.65) to (5.70) are known as the Kadanoff-Baym equations [7, 5, 8]. For a given ap- 
proximation of E[G] these equations can solved by time-propagation. In practice the equations 
are first solved for G M on the contour from to to to — i(3 parallel to the imaginary axis using the 
Kubo-Martin-Schwinger boundary conditions [9, 10, 11]. This amounts to a solution of the sta- 
tionary equilibrium problem. In fact since the Green function on the vertical part of the contour 
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corresponds to an equilibrium sitauation and only depends on r = T\ — r 2 the Eqs.(5.69) and (5.70) 
can be replaced by a single equation 

(-a r -/ l (l))G M (xi,x 2 ;r) = ^(r) + / M (xi,x 2 ,r) (5.76) 

where we defined 

G' M (x 1 ,x 2 ;r 1 - r 2 ) = G M (xm, x 2 r 2 ) (5.77) 

The functions such obtained can then be used as starting values for the real time propagation of 
the functions G^, and G"! . The initial conditions for these functions are then given by 

Gl(t„,T) = G M (0,r) (5.78) 

tf(r,to) = G m (t,0) (5.79) 

G<(t ,t ) = G M (0,0+) (5.80) 

G>(t ,t ) = G M (0+,0) (5.81) 



Exercise 



Derive the form of the Kadanoff-Baym equations given in Eqns.(5.65) to Eqns.(5.69) from the equa- 
tion of motion on the contour and check the initial conditions Eq.(5.78) to Eq.(5.81). 



The solution of the equations of motion is further simplified by the relations 

[G^(l,2)]* = -G^(2,l) (5.82) 
G>(xit,x 2 t) = -t<5(xi -x 2 ) + G < (x 1 i,x 2 i) (5.83) 

This means that in the time propagation one can restrict oneselves to solving G > (1,2) for t\ > t 2 
and G<(1,2) for t 2 < h. 



Exercise 

Prove relations (5.82) and (5.83) 



Let now discuss some physical properties that are contained in these equations of motion. As an 
example we write out Eq.(5.65) out in more detail. We then have 

(id tl -h(l))G < (l,2) = J dx 3 S HF (x 1 t 1 ,x 3 ti)G < (x 3 ti,2) 

+ / 1 d3[S > (l,3)-S < (l,3)]G < (3,2)- I * d3 £<(1, 3)[G>(3, 2) - G<(3, 2)] 

J to J to 

d3Gl(l,3)Sl"(3,2) (5.84) 
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The Hartree-Fock term on the right hand side of the equation is a potential that is spatially nonlocal, 
but local in time. Consequently it has no memory and does not lead to dissipation. The next to 
terms on the other hand involve time-integrations over all previous times. The complex self-energy 
kernels in these equations lead to dephasing and dissipative effects in the Green functions and lead 
to a decay of the Green function far from the time-diagonal. The last term on the right hand side 
is (apart from the Hartree-Fock term) the only term that remains for t\ = ti = to- It therefore 
describes the initial correlations in the system. 
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Chapter 6 



Conserving approximations 



6.1 Conservation laws and $-derivability 



There are several general exact relations known for systems in time-dependent external fields. These 
relations reduce in absence of the external field to the conservation laws for energy, momentum, 
angular momentum and particle number. The main question is now whether these conservation 
laws are also obeyed if we calculate the energy, momentum and angular momentum from the Green 
function obtained within a certain approximation. Approximations that do conserve them will be 
denoted as conserving approximations. The main question is then: how can we guarantee that a 
given approximation to the self-energy gives a conserving approximation for the Green function. 
Let us first discuss a number of conservation laws. If we know the Green function we can calculate 
the density and the current density from 

<n(l)> = -iG(l,l+) (6.1) 



{[£-¥ + H G(M#) } 



(6.2) 



r=i+ 

An important relation between the two quantities is provided by the continuity equation: 

S tl (n(l)) + V-<j(l)) = (6.3) 

This relation tells us that accumulation of charge in a certain region of space is related to current 
flow into that region. The is certainly an important relation that one wants to have satisfied in the 
analysis of charge transport processes. If we know the current density we can further calculate the 
total momentum and angular momentum expectation values in the system from the equations 

<P(*i)> = /^0(1)) (6-4) 
<L(ti)> = j dxin x (j(l)) (6.5) 



38 



For these two quantities the following relations should be satisfied 



a tl (P(t!)) = -y*dx 1 [(n(l))E(l) + (j(l)) xB(l 
^(Lfa)) = -ydxi[<n(l))n xE(l) + n x ((j(l)) xB(l)) 

where E and B are the electric and magnetic fields calculated from 

E(l) = Viu(l) - d tl A(l) 
B(l) = VixA(l) 



(6.6) 
(6.7) 



(6.8) 
(6.9) 



The equations (6.6) and (6.7) tell us that the change in momentum and angular momentum is equal 
to the total force and total torque on the system. In the absence of external fields these equations 
express momentum and angular momentum conservation. Since the right hand sides of Eq.(6.6) 
and Eq.(6.7) can also directly be calculated from the density and the current and therefore from 
the Green function, we may wonder whether they are satisfied for a given approximation to the 
Green function. Finally we will consider the case of energy conservation. Let E{t\) = (H(ti)) be 
the energy expectation value of the system, then we have 



d tl E(h) = - J dxi <j(l)> ■ E(l) 



(6.10) 



This equation tels us that the energy change of the system is equal to the work done on the system. 
Again we can ask whether this equation is satisfied for a given approximation to Green function. Let 
us first explain how the energy is calculated from the Green function. First of all, the one-particle 
energy is simply calculated from 



(/i(l)) = -i J dxi/i(l / )G , (l,l') i# 1+ = -« J dx 1 h(x 1 ,t 1 )G < (x 1 t 1 ,x' 1 t 1 ) 



x =xi 



(6.11) 



To calculate the expectation value of the intercation energy we simply use Eq.(5.1) to obtain 

id tl G<(l,2) = h(l)G<(l,2) + i J d3u;(l,3)<ft(2)ft( 3 )^tf(3)^tf(l)} (6.12) 

and therefore 
(id tl -h(l))G<(l,2) 



*2=tl 



/ 



C?X 3 l(;(xi,X3)(^(x2tl)^(x3ti)^//(x3tl)^//(xitl)) = 2i(W{t 1 )) 

(6.13) 



We the obtain for the total energy the expression 



E(t 1 ) = (h(t 1 )) + {W(t 1 )) 



dxi (id tl +/i(xi,ti))G < (xiii,x' 1 ti) 



(6.14) 



The question is now whether with this energy calculated form an approximate Green function and 
with the current density calculated from the same Green function relation (6.10) is satisfied. These 
questions were answered by Baym [12] in a famous paper. The main conclusion of the paper is as 
follows : If the self-energy £ is obtained from an underlying functional $ [G] , according to 



£(1,2) 



<5$ 
<5G(2,1)' 



(6.15) 
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+ 




Figure 6.1: Some of the low-order diagrams, and some of the self-energy diagrams obtained from 
£ = 5&/5G. The prefactor of a ^-diagram is ns/2n where ns is the number of topologically 
different S-diagrams that can be generated from it and n is the number of interaction lines. 

then solving the Kadanoff-Baym equations with this approximation to the self-energy will lead to 
a Green function G that satisfies exactly the relations (6.3), (6.6), (6.7) and (6.10). The question 
is now how the functional can be constructed. Such a functional can be constructed, as first 
shown by Luttinger and Ward [13], by summing over irreducible self-energy diagrams closed with 
an additional Green function line and multiplied by appropriate numerical factors, 



The term n indicates the number of interaction lines and k labels S-diagrams. The trace tr 
indicates an integration over all variables (in contrast to the trace Tr that denotes a summation 
over a complete set of states in the Hilbert space). Some of the low-order diagrams are shown in 
Fig. 6.1, together with some of the corresponding self-energy diagrams. 



6.2 Approximate conserving schemes 

6.2.1 Time-dependent Hartree-Fock approximation 

We now consider some approximate conserving schemes. One of the simplest approximations we can 
take is the Hartree-Fock approximation for E. For this approximation the ^-diagrams correspond 
to the first two in Fig. 6.1 and the self-energy is explicitly given by Eq.(5.52) or in more detail in 




(6.16) 
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Eqs.(5.63) and (5.64) In that case the equations of motion for G^ are given by 

(id tl - h(l))G$(l,2) = ydx 3 S HF (l,X3ti)G^(x 3 ti,2) (6.17) 

{-idt a - h(2))G$(l,2) = ydx 3 G^(l,X3t 2 )S HF (x 3 t 2 ,2) (6.18) 

To solve these equations we first introduce the following orbitals 

(id tl - h(l))tpi(l) = y'dx 3 S HF (l,x 3 t 1 )( / 9 J (x 3 t 1 ) (6.19) 

(-id t2 - h(2))tpi(2) = ydx 3 ^(x 3 t 2 )S HF (x 3 t 2 ,2) (6.20) 

From the basic property (5.82) and the explicit form of S HF of Eq.(5.64) it follows that T, HF is 
hermitian and therefore (pi = (p* (this is not true anymore on the vertical part of the contour) . We 
now make the following Ansatz for G < : 

G<(l,2) = i5] W (%;(2) (6.21) 
j 

where rij are constants to be determined. It is clear that this Ansatz for G < satifies Eqs.(6.17) and 
(6.18) but we still need to check the boundary conditions. In order to satisfy Eq.(5.83) we make 
for G > the Ansatz 

G>(1,2) = -^(1-^)^(1)^(2) (6.22) 

3 

such that 

G>(1,2) -G<(1,2) =-iViP;(%J(2) =-i(5(x 1 -x 2 ) (6.23) 



ti=t 2 

3 



ti=t 2 



We must further specify the initial conditions. This can be done with help of Eqs.(5.80) and (5.81). 
We therefore first have to solve for the Green function G M on the imaginary part of the contour 
where t = —it as given in Eq.(4.13). 

G M (xiri,x 2 r 2 ) = 0(n -r 2 )G > (x 1 ,-ir 1 ,x 2 ,-ir 2 ) +^(r 2 -ri)G < (xi,-ir 1 ,x 2 ,-ir 2 ) (6.24) 

If we define Xi( x i r i) = ¥>i( x ii — * r i) an d Xi( x i T i) = ¥>* ( x i> ~~ * T i) then we have 

G > (x i , -iri,x 2 , -ir 2 ) = -i^(l - n j )xj(xiri)xj(x 2 r 2 ) (6.25) 

j 

G < (x i , -iri,x 2 , -it 2 ) = i^n j Xi(xiri)xj(x 2 r 2 ) (6.26) 

j 

where Xi an d X« satisfy the equations 

(-d Tl - /i(xi))xi(xin) = y'dx 3 S HF (x 1 r 1 ,x 3 r 1 )^(x 3 T 1 ) (6.27) 
(<9 T2 - /i(x 2 ))xi(x 2 T 2 ) = y , dx 3 Xi(x 3 r 2 )S HF (x 3 T 2 ,x 2 r 2 ) (6.28) 
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With these equations it is readily seen that the equations of motion for G M of (5.70) and (5.76) 
are indeed satisfied. We therefore only have to specify the boundary condition. On the vertical 
contour the system is in equilibrium and therefore E HF (xir, x 2 t) does not depend on r. For the 
functions Xi an d Xi we can therefore make the Ansatz 

Xi(xiri) = ^(xi)e- £m (6.29) 

Xi(xin) = ^(xije 6 * 71 (6.30) 

If we insert these forms in Eq.(6.27) and (6.28) we obtain 

(€i ~ /»(xi))^i(xi) = y ( iX3S HF (x 1 T 1 ,X3T 1 )^(x3) (6.31) 

- / i (x 2 ))V*(x 2 r 2 ) = y , dx3V*(x3)S HF (x 3 r 2 ,x 2 r 2 ) (6.32) 

and we obtain 

G > (x 1 ,-ir 1 ,x 2 ,-ir 2 ) = - n,)e-^ (ri - r2) Vi(xi)^(x 2 ) (6.33) 

j 

G < (x 1 ,-ir 1 ,x 2 ,-ir 2 ) = i£ ^-^-^(xi^x,) (6.34) 



From these equations and Eq.(5.64) we see indeed that on the imaginary part of the contour the 
self-energy is time-independent and given by 

S HF (xi,x 2 ) = -u>(xi,x 2 ) ^njY>i(xi)V>*(x 2 ) +<5(xi - x 2 ) J dx 3 io(xi, x 3 ) ^ nj^(x 3 )V>*(x 3 ) 

j j 

(6.35) 

Let us now determine the coefficients rtj from the Kubo-Martin-Schwinger boundary conditions. 
From condition Eq.(4.18) and the explict forms (6.33) and (6.34) we see that we must have 

G < (xi,0,x 2 , -it 2 ) = i^n i e £lT2 ^ j (xi)V'*(x 2 ) 

j 

= - % -)e" ei(/3_T2) V' J -(xi)^(x 2 ) = -G > (x 1 ,-i/3,x 2 ,r 2 ) (6.36) 

We thus obtain the relation 

n, = (l-n,)e-^ (6.37) 

this equation is easily solve to give 

1 

nj = -s 6.38 

eP*i +1 v ' 

which is the famous Fermi-Dirac distribution. We have now also completely determined the initial 
conditions of the orbitals ip, we have <£j(xO) = ^j(x) where ipi solves Eq.(6.31) with self-energy 
(6.35). Now that we have completely determined we can also explicitly evaluate the self-energy 
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of Eq.(5.64). The results are summarized below: 

G<(1,2) = ^n^(l)^(2) (6.39) 
j 

G>(1,2) = -^(1-^)^(1)^(2) (6.40) 

(i^ - = y^x 3 S HF (l,X3ti)^(x3ti) (6.41) 

S HF (xiti,x 2 ti) = -w(xi, x 2 ) ^ njy>i(xiti)y*(x 2 ti) 

i 

+ <5(xi-x 2 ) / rfx 3 w(xi,X3)y^njyj(x3ti)y*(x 3 ti) (6.42) 
^ j 

These equations are just the time-dependent Hartree-Fock (TDHF) equations. If we choose the 
chemical potential fi between the highest occupied and lowest unoccupied level of the stationar 
Hartree-Fock equations and take the zero temperature limit we find that nj becomes equal to one 
for the occupied states and zero for the unoccupied states. 

It is now also interesting to see what the time-dependent Hartree-Fock approximation will give 
for the density response function. From Eq.(5.48) we see that the first order change in the Green 
function due to a change in the external field is given by 

5G(l,l') = J d3d4G(l,3)G(4,l')T(34;2)6v(2) (6.43) 

and therefore the density reponse is given by 

Sn(l) = -iG(l,l + ) = ~ij d3d4G(l,3)G(4,l)T(34;2)6v(2) = J d2 x (l,2)6v(2) (6.44) 

and the density response function has the expression 

X(l,2) = ~ij d3d4G(l,3)G(4,l)r(34;2) (6.45) 

The density response function corresponding to TDHF is then obtained by inserting the vertex 
function V that belongs to TDHF. This function is then the solution of Eq.(5.61) where we need to 
insert (5S HF (1, 2)/<5G(4, 5) for the kernel. We have 

SG( l = 2)6(1, 4)5(2, 5) - iS(l, 2)6(4, 5+)w{l, 4) (6.46) 

and therefore from Eq.(5.61) we obtain 

r(12;3) = 5(1,2)5(1,3) + iw (1,2) J d6d7 G (1,6) G (7, 2)T (67; 3) 

-i6(l,2) J d5d6d7w(l,5)G(5,6)G(7,5)T(67;3) (6.47) 
This gives indeed the well-known 'ladder and bubbles' series for the TDHF response function. 
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6.2.2 Second Born approximation 

The level conserving approximation contians the second order diagrams, corresponding to the first 
four <&-diagrams of Fig. 6.1 This corresponds to the so-called second Born approximation. Let the 
second order diagrams be denoted by X B then there explicit form is given by the last two terms of 
Eq.(5.56), i.e. 

£(1,2) = £ HF (1,2) + £ B (1,2) (6.48) 
£ B (1,2) = i 2 j d3d4G(l,3)u>(l+,4)G(3,4)G(4,2)w(3 + ,2) 

- i 2 J d3d4G(l,2)u/(l + ,3)w(2,4)G(4,3)G(3,4 + ) (6.49) 

Because of the time delta functions contained in w(1.2) we can write S B as 

S B (1,2) = i 2 J dx3dx4G(xiti,X3t2)w(xi,X4)G(x3t2,X4ti)G(x4ti,X2t 2 )w(x3,X2) 

- i 2 J dx 3 dx 4 G(xiti, x 2 t2)w(xi, x 3 )ti;(x2, x 4 )G(x 4 t 2 , x 3 ti)G(x 3 ti, x 4 t 2 ) (6.50) 

We see that these expressions do not involve and time integrations. Therefore it is straightforward 
to find the various S^-components of this expression where for x we have >, <, [~, ] , M. We have 

S B '^(1,2) = i 2 J dX3dx 4 G^(xiti,X3t 2 )-a)(xi,X4)G^(x3t 2 ,X4ti)G^(x4ti,X2t2)^(x3,X2) 

- i 2 J dx3dx 4 G^(xiti,X2t2)w(xi,X 3 )?i;(x 2 ,X4)G^(x4t2,X3ti)G^(x3ti,X4t2) (6.51) 

£ B 'l (1,2) = i 2 J dx 3 <ix 4 G" 1 (xiti, x 3 r 2 )tt;(xi, x 4 )G r (x 3 T2, x 4 ii)G _l (x 4 ii, x 2 T 2 )iy(x 3 , x 2 ) 

- i 2 j dx 3 dx4G^(xiti,X2r2)ti)(xi,x 3 )ti;(x2,X4)G r (x4r2,X3ti)G^(x3ti,X4r 2 ) (6.52) 
S B 'T(1, 2) = j2 j dx 3 dx A G^xin, x 3 i 2 )w(xi, x 4 )Gl (x 3 t 2 , x 4 ri)G r (x 4 ri, x 2 t 2 )w(x3, x 2 ) 

- i 2 j dx 3 dx 4 G r (xiTi,x 2 t2)^(xi,X3)ti;(x2,X4)G" l (x4t2,X3ri)G r (x3ri,X4t2) (6.53) 
S B - M (1,2) = i 2 J dx 3 dx4G M (xiri,X3r 2 )u;(xi,x 4 )G M (x3r2,X4r 1 )G M (x 4 ri,x 2 r2)u;(x3,X2) 

- i 2 j dx3dx 4 G M (xiri,X2r2)w(xi,X3)ti;(x2,X4)G M (x4r2,X3Ti)G M (x3ri,x 4 r2) (6.54) 
With these expressions and the Kadanoff-Baym equations the problem is completely defined. 



Exercise 

Suppose we calculate the first order density change by propagation of the Kadanoff-Baym equations 
within the second Born approximation. What would be the diagrammatic structure of the reponse 
function that we obtain in this way ? 
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6.3 Applications 



In this section we will illustrate the methods discussed sofar with some results. When one does 
actual calculations the Green functions are often expressed in a basis, i.e. one writes 

G(xiii,x 2 t2) = y^yi(xi)y*(x 2 )G^(ti,t2) (6.55) 



where <pi represents a suitable chosen basis such as Hartree-Fock molecular orbitals. The coefficients 
Gij are in fact the Green 
with respect to this basis 



Gij are in fact the Green functions with respect the annihilation and creation operators a, and a\ 



GijitxM) = -i(Tc[hH{h)a> hH (t 2 )]) (6.56) 
This means that an equation like (5.84) attains the form 

(id tl -h(t 1 ))G < (t 1 ,t 2 ) = V UF (t 1 ) ■G < (t 1 , t 2 ) 

+ I 1 dt 3 [£>(*!, t 3 ) - *s)] • G<(t 3 ,t 2 ) - [ 2 dt 3 S<(t u t 3 ) ■ [G>(t 3 ,t 2 ) - G<(t 3 ,t 2 )} 

J to J to 

r/3 

- i / dT 3 G\t U T 3 )^(T 3 ,t 2 ) (6.57) 
Jo 

where all be bold-faced symbols now represent matrices and the symbol " • " a matrix product. 
As a first application we consider a laser-excited quantum well of GaAs [14] where we only take the 
valence and conduction bands as basis functions. For this case the Kadanoff-Baym equations can 
be solved within the second Born approximation. The laser field is taken to be a 50 fs laser pulse 
with its maximum peak intensity at t = 0. In Figs. (6.2) and (6.3) the imaginary parts of Gf c and 
Gy V are displayed (both taken at the top of the valence band, i.e. in the state k = 0). As follows 
from Eq.(5.82) these quantities are indeed symmetric in t\ and t 2 . In general one has 

ImGf j (t 1 ,t 2 ) = lmGf i (t 2 ,t 1 ) (6.58) 
Let us see what this quantity represents. Now since 

Gf j {t 1 M) = i{a\ H {t 2 )a i ,H{ti)\) (6-59) 

it follows that 

ImGuiUM) = (dl H (ti)di,H(ti)} = riiih) (6.60) 

can be interpreted as the occupation number of the state i. This is indeed consistent with figures 6.2 
and 6.3. One sees that due to the laser pulse the condition band gets occupied (the time diagonal 
in Fig. 6.2) whereas the valence band gets depleted (the time diagonal in Fig. 6.3). One further sees 
that the functions decay away from the time-diagonal. This is due to to the memory terms in the 
Kadanoff-Baym equations, i.e. the second and the third term on the right hand side of Eq.(6.57). 
The decay would be absent if we would only take into account the Hartree-Fock term. 
As a next application we solved the Kadanoff-Baym equations within the second Born approxima- 
tion for a hydrogen molecule. The results are very recent and therefore we only plot in Fig. 6.4 
the function Im G^ gCFg (ti, t 2 ) for a field free propagation without a laser pulse excitation. We see 
a similar picture as for the two-band GaAs calculations. The function oscillates away from the 
diagonal with a frequency that is close to the ionization energy (as it should). 
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Figure 6.2: Conduction band part Im Gf c {t\, t^). A 50/s laser pulse is applied to GaAs which has 
it peak strenght at t = O.The time-diagonal gives the occupation number of the conduction band. 
The decay away from the time diagonal is due to electron correlations. 
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Figure 6.3: Valence band part Im £2)- A 50/s laser pulse is applied to GaAs which has it 
peak strenght at t = O.The time-diagonal gives the occupation number of the valence band, it is 
clearly seen that the laser excitation leads to depletion of the valence band. The decay away from 
the time diagonal is due to electron correlations. 
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Figure 6.4: Im G agCTg (t\ , tz) for the H2 molecule obtained from field free propagation of the Kadanoff- 
Baym equations within the second Born approximation 
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Chapter 7 



Outlook 



In these notes we have given an introduction to the nonequilibrium Green function method, starting 
from basic ideas based on the time contour to the final derivation of the Kadanoff-Baym equations 
which are the ones that need to be solved in practical applications. Due to the limited space 
necessarily many things have been left out. Many things can still be said about the properties 
of the spectral functions that can be calculated from the nonequilibrium Green function method 
and which are related to other exact identities known as the Ward identities. Also the topic of 
more general initial states than the equilibrium state has been left out. The properties of the 
^-functional have only been briefly touched and other conserving schemes such as the GW-method 
and the T-matrix approximation have not been discussed at all. 

Nevertheless we believe that the background presented here will provide the reader with enough 
information to study the wide literature on transport throught single molecules, as well as many 
other applications of nonequilibrium Green function theory. Especially with the developments in 
single molecule conduction the interest in the nonequilibrium Green function method has been 
growing steadily and is likely to do so in the near future. 
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